arXiv: 1509.04548vl [quant-ph] 15 Sep 2015 


Correlations of photon trajectories in the 
problem of light scintillations 

R. A. Baskov and O. O. Chumak* 


Institute of Physics of the National Academy of Sciences 
pr. Nauki 46, Kyiv-28, MSP 03028 Ukraine 


September 16, 2015 


Abstract 

A distribution function approach is applied to describe the dynam¬ 
ics of the laser beam in the Earth atmosphere. Using a formal solution 
of the kinetic equation for the distribution function, we have devel¬ 
oped an iterative scheme for calculation of the scintillation index (cr^). 
The problem reduces to obtaining the photon trajectories and their 
correlations. Bringing together theoretical calculations and many-fold 
computer integrations, the value of is obtained. It is shown that 
a considerable growth of in the range of a moderate turbulence is 
due to the correlations of different trajectories. The criteria of appli¬ 
cability of our approach for both the coherent and partially coherent 
light are derived. 


1 Introduction 

Basic principles of the radiation transfer theory were formulated in the sem¬ 
inal paper of Schuster [1] as early as the beginning of the 20th century. The 
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paper pQ was devoted to the light propagation in a foggy atmosphere. Since 
then the Schuster approach has obtained many applications in such impor¬ 
tant helds as astronomy, laser communication and radar systems, remote 
sensing etc. 

The range and performance of light communication systems are limited 
signihcantly by an unfavorable influence of local fluctuations of the refractive 
index of the Earth atmosphere. On the other hand, high sensitivity of the 
photon trajectories to the fluctuations can be used for the atmosphere diag¬ 
nostics [2] . The key issue about the index-of-refraction structure constant 
is that it cannot be reliably computed from Erst principles. The refractive- 
index fluctuations arise from temperature inhomogeneities of the air. The 
inhomogeneities cause turbulent eddies which give rise to a random distribu¬ 
tion of the air density([3]-|5]). This results in random spatial variations of 
the refractive index. 

The turbulent eddies are described by a wide range of characteristic 
lengths of inhomogeneities. These lengths cover the interval from few mil¬ 
limeters (the inner radius, Iq) to hundred meters (the outer radius, Lq). 
Therefore, various types of beam scattering are observed. The scattering 
by large-size eddies results in random redirections of the beam as a whole. 
This process is known in the literature as a ’’wandering” or ’’dancing” of the 
beam 0 . 0 . On the other hand, the scattering by small-size eddies causes 
spreading of the beam. For a long-distance propagation or a strong turbu¬ 
lence, the beam radius becomes greater than the characteristic sizes of the 
inhomogeneities. In this case the probability of the beam to be redirected 
becomes small and the relative value of the wandering radius decreases |8]. 

The beam wandering and broadening can be considered as the specihc 
manifestations of a more general phenomenon, namely the intensity fluctua¬ 
tions (i.e. scintillations) caused by the atmosphere turbulence. The scintil¬ 
lations have a tendency of saturating for a long-distance propagation 0.110 
(the regime of a strong turbulence). This is because in the course of propaga¬ 
tion the radiation acquires the properties of the Gaussian statistics when the 
signal-to-noise ratio (SNR) tends to unity. The asymptotic behavior of the 
scintillation index, —)■ 1, was explained in Refs. nn-na. Moreover, it was 
shown quite generally that this property stays unchanged for any refractive 
index distribution, provided the response time of the recording instrument is 
short compared with the source coherence time. This result was conhrmed 
analytically in [IT] . 

At the same time, calculations, performed by different methods in [T5] 
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and [H], show a possibility of significant snppression of the scintillations. 
To this end partially coherent laser beams with the coherence time shorter 
than the detector integration time (a slow detector) can be nsed. The case 
of a partial coherence was also stndied in ini .|18]. Recent theoretical and 
experimental developments on propagation of partially coherent beams in a 
tnrbnlent atmosphere were discnssed in [12]. 

There are several analytical approaches explaining behavior of the scin¬ 
tillation index in the case of strong tnrbulence [1611201121]. Their analysis is 
based on the physical pictnre where fonr waves, forming the second moment 
of the intensity, conserve only pair correlations in conrse of long-distance 
propagation. Two different pairs of the photon trajectories contribnte into 
the sqnare of the photon density at the detector. Dashen nsed the Feynman 
path integrals to prove that in a convincing manner [20] . 

The recent interest to beam propagation was awakened by the develop¬ 
ment of qnantnm commnnication in the free atmosphere [22], [23] . Detailed 
stndies of the effect of the tnrbnlence-indnced losses on the qnantnm state of 
the light in the conrse of satellite-mediated commnnication and for realiza¬ 
tion of the entanglement transfer in the atmosphere were reported in Refs. 
[21] and [25]. 

The formalism of the photon distribntion fnnction (the photon density in 
the coordinate-momentnm space [26]), is also applicable to the problem of 
scintillations [SI [m EH ES]. The mentioned papers are based on a physical 
pictnre, which is similar to the described above. The method of photon distri¬ 
bntion fnnction is nsed for description of both the classical and the qnantnm 
light inclnding propagation of single-photon pnlses (see, for example. Refs. 
[271 I29l [30]). Solntion of the kinetic eqnation for the operator of photon 
density is based on the method of characteristics. The assnmption of weak 
distnrbances of photon momenta by the atmosphere (the paraxial approxi¬ 
mation) rednces the problems of scintillations to the problem of obtaining 
photon trajectories and their correlations. A slowly varying flnctnating force, 
deflecting photon trajectories from straight lines, describes the effect of the 
atmospheric eddies. 

In this work we stndy the scintillation index for moderate and strong 
tnrbnlences, when correlation of trajectories of only two photons is reqnired. 
Accnracy of the calculations depends on the accuracy of obtaining the tra¬ 
jectories. Using high-order iterations and bringing together analytical and 
numerical procedures, we calculate the scintillation index. Our main interest 
is to analyze the range of moderate turbulence strengths where previous the- 
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ories do not ensure a reliable description. Comparison of the obtained results 
with those represented in [TB] helps indicate the range of turbulence where a 
simplihed approach should be corrected by high-order iterations. Also, our 
studies describe more realistically the effect of partial coherence. 

2 Photon distribution function approach 

The photon distribution function is dehned by analogy with distribution 
functions in solid state physics. In particular, it is similar to the phonon 
distribution function. Both of them are dehned as [26].[3T] 

/(r, q, ^ (1) 

k 

where and 6q are the bosonic creation and annihilation operators of pho¬ 
tons or phonons with the momenta q, and V = L^LyL^ is the normalizing 
volume. Polarization of the corresponding modes is not specihed in ([1]). 
In the paraxial approximation, assumed here, the initial polarization of the 
beam remains almost unchanged even for a long-distance propagation (see, 
for example. Ref. [32]). 

The operator /(r,q, t)) describes the photon (phonon) density in the 
phase (r,q) space. Usually, the characteristic sizes of spatial inhomogeneities 
of the radiation held are much greater than the wave-length. In this case 
the sum in Eq. ([2]) can be restricted by small k. Here and in what follows 
we consider that k < ko qo, where go is the wave vector corresponding to 
the central frequency of the radiation, uq = cqq. At the same time /cq should 
be taken sufficiently large to provide a required accuracy of the beam prohle 
description. 

The evolution of the Heisenberg operator /(r, q, f) is determined by the 
commutator 

^t/(r,q,t) = ^[/(r,qU),h^], (2) 

where 

H = '^ hiUqfot 6q - ^ hWqUkfeijfeq+k, (3) 

q q,k 

is the Hamiltonian of photons in a medium with a huctuating refractive index 
n(r) (rik is its Fourie transform), = hcq and Cq = |^ are the vacuum 
values of the photon energy and velocity, respectively. 
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Assuming the characteristic values of the photon momentum to be much 
greater than the wave vectors of turbulence, the kinetic equation for the 
photon distribution function can be written as 


+ Cqar+ F(r)aq}/(r,q,t) = 0, (4) 

where F(r) = uodruir) is the random force originating from the atmospheric 
turbulence. The general solution of Eq. (jl]) is given by 


/(r,q,t) 




, Mt') 




(5) 


where the function 0(r, q) is the ’’initial” value of /(r,q, f), i.e. 

0(r, q) = -^ Vk/ 2 ) \t=o = q)- (6) 

k k 


The derivatives 


dr{t') 

at' 


and 


dp(t') 

at' 


should satisfy the equations 


dr{t') 

dt' 

dt' 


c[q(f')] 

F[r(T)], 


(7) 


completed with the boundary conditions r(f') = r and q(T) = q for t = t'. 
As we see, Eqs. ([7]) coincide with the classical (the Newton) equations of 
motion of a point particle moving with the velocity Cq and affected by an 
external force F(r). Formal solutions of Eqs. o can be written as 


q(T) = q + ^ df"F[r(q,f")] 


and 


r(T) = r - Cq(f -t')- — [ - t')F[r{t")], 

Qo Jt 

Eqs. ([H]) and dH]) allow us to rewrite the expression (jS]) as 

/(r,q,t) = r - Cqt + — /" dt't'F[r{q,t%ci- [ dt'F[r{q,t' 
[ Qo Jo Jo 


)] • 


( 8 ) 

(9) 


( 10 ) 
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If F(r) is a known function, an approximate value for /(r, q, t) can be ob¬ 
tained by inserting the term r(q, t') r — Cq(t — t') into Eq. fllUl) . In this 
case the argument of the fluctuating force F[r(q, t')] is replaced by a straight 
line, that is correct only in the absence of the turbulence. Improvement of 
the theory can be achieved if the argument of F accounts for the turbulence. 

It follows from Eq. flTU]) that statistical properties of the radiation depend 
not only on the turbulence but also on the initial distribution function 0(r, q). 
This function is determined by the source held. Its explicit form is determined 
in the course of ’’sewing” of the near-aperture and the atmospheric helds [16] 
given by the amplitudes &q(&q). We consider the light propagation in the 
z-direction. The source held is assumed to be described by the Gaussian 
function, ‘h(r) = (2/7r)^/^Tg->'_LAo. Then the propagating amplitudes are 
given by 

= 0 ) = ( 11 ) 

where b is the near-aperture amplitude of the laser held, index (x) means the 
perpendicular to the ^-axis components, and S = L^Ly. 

We will take into account the ehect of the phase dihuser by multiplying 
the distribution <I)(r) by the phase factor where the quantity a is a 

random variable. In this case Eq. flTT]) should be modihed by substituting 
in its right-hand side q_L -|- a = qa for q^. Such a simple modeling of the 
phase dihuser is justihed if (i) the detection time is much longer than the 
characteristic time of the variation of a (slow detector) and (ii) there is a 
large root mean square of the phase huctuations. (More detailed analysis is 
presented in [2H]-) This case corresponds to the Gaussian distribution of a: 


P{,0,x,y) 


A 

2 ^ 




( 12 ) 


with a covariance {ai „) = X ^ and the transverse correlation function of the 
outgoing held (at t = 0) is given by 

(F(r^)F(r^ + A))^ = (13) 


Here Eq = E{r± = 0,A = 0,t = 0) and the notation (...)» means averag¬ 
ing over distribution P{ax,y)- The radiation, whose correlation properties 
are described by function ffT3|l . is referred to as the Gaussian Shell-model 
held. The parameter A in the exponential factor describes the decrease of 
the transverse correlation length. It can also be said that this parameter 
generates a new characteristic length, 1/ri, in the momentum distribution 
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(i.e., in the q-domain). This is seen from the explicit term for 0(k, q) which 
after averaging over the fluctuations of a reduces to 


{4){k, q))a 


b^b 


2 -(?■ ^ 

2 ''x 8 ^ 


(14) 


where rj = Tq + 2rgA ^ j , and variables and kz are omitted. 

It is seen from Eq. fflTD that qj_ is distributed in the range of the order 
of \/2lri that is greater than the one for coherent beam. In contrast, the 
characteristic value of k depends only on the initial size of the beam ( k ~ 
x/8/ro). 

In the course of light propagation, the diffraction phenomena and scatter¬ 
ing by atmospheric inhomogeneities broaden the beam resulting in decrease 
of k. At the same time, the value of q increases with the distance. This is 
because of the Brownian-like motion of photons in the q^-domain (see Ref. 
[T6]). Such a simple physical picture, elucidating evolution of the beam ge¬ 
ometry, is, however, not applicable to the description of scintillations. The 
phenomenon of scintillations is more complecated and can be described in 
terms of spatio-temporal correlations of four waves. 


3 Scintillation index 


The photon distribution function is used here to obtain the scintillation index 
(T^. The dehnition of cr^ is given by 


, {P{v)) - (/(r))^ 

{m? 


(15) 


The photon density I{r,t) is expressed in terms of the distribution function 
as 


q 


27r- 


b'^brQ 

'W 


^-ik[r-c(q)t+^ f* dt't'F(r(q,t'))]-Ql^-k^^ 

q,k 


t (16) 

where Qa = Q-(-a=q-|-a — dTF[r(q, t'). The summation is taking over 
qx and components, while and kz are considered to be hxed: = go 

and kz = 0. The exponential term originates from the solution ffTOj) of the 
kinetic equation (jl]). 
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To obtain (/(r,t)), three independent averagings are reqnired. One of 
them concerns the sonrce variables. In the case of a coherent state of the 
source, |/3), we have {Wh) = |/9p. The second averaging over a random phase 
of the diffuser should be carried out as explained by Eq. ffTT|) . The third 
averaging deals with the fluctuating force F. These three actions can be 
performed independently that facilitates the analysis. Also, the calculations 
are simplihed if we use the identity 


,-QV2/2 ^ 


dp 

27rrg 


^ipQ-p2/2r2 


(17) 


Because of Eq. (HU, the term in the exponent of Eq. (1T6|1 reduces to the lin¬ 
ear in F form. Then, considering F as a random Gaussian variable, the value 
of (/(r, t)) can be easily obtained in a manner, explained in Ref. [16]. To 
calculate (/(r, t)), an explicit form of the refractive-index correlation func¬ 
tion, (n(r)n(r')), is required. In a statistically homogeneous atmosphere it 
can be written as 


(n(r)n(r')) = / dge '’'V(g) 


(18) 


A widely used the von Karman approximation for the spectrum, '0(g)) is 
given by 


0(g) = 0.033^2 


exp[-(^/o/27r)^] 


g| = d, 


(19) 


where the vector g is dehned in the three dimensional domain. 

The ’’source” part of (/^(r)), given by (b^bb^b), is approximately equal to 
{b^b^bb) = 1/50, when the condition |/50 >> |00 is satished. This inequality 
implies that the initial laser radiation is in a multiphoton coherent state. The 
averaging over independent random quantities a and a' can be used instead 
of the time averaging of the diffuser state. Then we have 


(/^(r,t)) = 


27r/52r2 




q,k,q',k' 


—ikfr—c, 




X + e-[(Q«+l0+(Q;-^0+(Qa'-|d+(QL,+^0]^}). ( 20 ) 

There are two terms in the braces of Eq. fl20|) . They appear only if the 
initial four-wave correlation reduces to the pair correlation [16] . Such a mod- 
ihcation of the statistical properties of the radiation occurs when the waves 
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.}S'l2 , /o' 








propagate for a long time which is sufficient for randomization of the trans¬ 
verse photon momentum. A more general case, which includes the regime of 
fast detection, was analyzed in Ref. [28] . 

The averaging of Eq. fl20|) over a and a' results in 


2Tif5^ri 


VS 


E 


^-i{k[r-c(q)i]+k'[r-c(q')t] + ^/o*rft't'[kF(r{q,t'))+k'F(r(q',t'))]} 


1 , 1 ,/ 


X |g-(Q^+Q''W/2-(fc2+fc'hrg/8_^g-[(Q-Q')2+(k+k')2/4]r=/4-[(Q+Q')2+(k-k')2/4]rf/4j^_ 

( 21 ) 

In the absence of a phase diffuser, ro = ri, the summands in the last braces 
contribute equally into ([21]). 

2 /2 

Similarly to Eq. di. the factor e ^ ^ in fl?I|) can be expressed in 

the integral form as 


g-{Q2+Q'2) 


ipQ+VQ'-(p"+p'h/2r-? 

(27rr^)2 


( 22 ) 


As we see, the exponent in the left-hand side is represented as a linear form 
of the force F. A similar transform is applicable to the second term in the 
last braces of fl2T]) . As a result, the fluctuating force enters the right-hand 
side of fl2T]) only via the common multiplier, M, given by 


][f — g-*/o <^*'{(P+kt'c/ 9 o)F[r(q,t')]+{p'+k't'c/go)F[r(q',t')]} 


(23) 


Obtaining of the average value of P reduces to averaging of M with many-fold 
integration. Assuming the exponent in (|2^ as a Gaussian random variable, 
we can write 


Two types of the correlation functions determine (M): 


— (}~^Vpp+‘^4>ppi+4>pip/) 


(24) 


t pt 


Oppi — 


jpp — 



'0 ^0 
pt pt 



dt'dt''{p + kt'c/qo) ■ (F[r(q, T)]F[r(q', t")]) ■ (p' -h k't''c/qo)], 

(25) 

dt'dt”{p + kt'c/ go) ■ (F[r(q, f')]F[r(q, t")]) ■ (p -F kfc/ go)], (26) 


0 JO 
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where symbols P and P' denote sets of three vector variables P = {q, p, k} 
and P' = {q',p',k'}. The correlation fnnctions of the forces along different 
(q 7^ q') and coinsiding (q = q') trajectories enter Eqs. fl25|) and fl26|) . 
respectively. The former can be rewritten as 

(F„[r(q,f')]^/3[r(q',t")]) = (^a[r(q,f') - r(q', f")]T>[0]), (27) 

where the notations a and /3 stand for the x and y - components. The 
expression for fl26l) follows from Eq. fl27|l by setting q = q'. 

The right-hand side of Eq. fl27)) is assnmed to be a fnnction of the coor¬ 
dinate difference, r(q, T) — r(q',f"). It is so if the atmosphere is statistically 
homogeneous. In the course of averaging, dependence of the coordinate dif¬ 
ference on the fluctuating force should be also taken into account. This 
dependence is given by the relation 

r(q,T)-r(q',f") = (e^c-fCq/)(T-t")-Cq-q/(t-f')H-/ dti(T-ti)F[r(q', fi)] 

% Jt' 

+ —f dti{t'-ti){F[r{q,ti)] (28) 

Qo Jt 

which follows from Eq. ([9]). The distance |r(q, T) — r(q',f")| should be of 
the order or less than the outer radius, Lq, of the turbulence. Taking into 
account that c » |cq_q'|, |cq/|, we infer that \t' — t''\ < Lq/c. This means 
that in the right-hand side of Eq. fl28|l Cq/ in the first term and the third 
term, which is proportional to {f — t")^, can be omitted. Then Eq. fl25]l 
reduces to 

r(q, t') - r(q', t") = e^c{t' - f) - Cq_q/(t - t') 

+ —f dti{t'-ti){F[r{q,ti)]-F[r{q^,ti)]}. (29) 

Qo Jt 

The last two terms in Eq. fl2^ describe the displacement of two photons 
from each other because of the difference of their initial velocities. The term 
—Cq_q/(f — f) describes the divergence of two straight-line trajectories. The 
last term accounts for the different actions of the atmosphere on the particles 
moving in different spatial regions. 

Obtaining of the average values in Eq. fl27)) . which depend on the wave- 
vectors r(q, f') and r{q',t"), seems to be challenging because of the presence 
of the fluctuating force in r(q, f') and r{q',t"). Nevertheless the analysis 
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simplifies if we neglect the correlations between the forces Fq, or and the 
forces entering r(q, t') or r(q', t”). This simplihcation can be justihed by the 
following reasonings. The explicit value of the a-force is given by 

F„[r(q,T)] = Fo[r - Cq(t - T)-/ dti(ti - t')F[r{q,ti)]] 

Qo Jt 

c 

= Fa[r±-+ cej,' -/ - t')F[r(q, F)]], (30) 

Qo Jt 

where the relation z = ct is used. 

If the correlation exists, the distance |r(q, ti) — r(q, f')| can be estimated 
by the value c(ti — t') < Lq. In this case, the integral in Eq. (1301) is propor¬ 
tional to {Lq/cY'. Hence, the correlation between FQ,[r(q,f')] and F[r(q,fi)] 
can be neglected. This approximation implies the physical picture where the 
variation of the photon momentum on the correlation length, Lq, is much 
smaller than qq. Therefore, the averaging (FaFp) can be performed in two 
steps. Firstly, we obtain {FaFp) considering the arguments of Fa and Fp 
to be fixed. After that, the averaging of the forces, entering the arguments, 
should be performed. For example, the term (|25ll is expressed as 


t ft 


bppi — cUq 



0 JO 


dtJdf / dg?/^(^)g-(p+kT-) 

J Qo Qo 


(31) 

where the hrst-step averaging results in appearance of the spectral density 
'ip{g). The second-step averaging is shown in flHTD by the angle brackets. To 
simplify the derivation of the authors of [IB] represented the average of 
the exponential function in Eq. fl3T|) as a product. 


(e 


-ig[r{q,t')-r{q',i' 


(e 


-igr(q,t')^^pigr{q',r 


">). 


(32) 


neglecting the correlation of the photon displacements r(q, t') and r(q',f"). 
Further analysis explains how this correlation can be accounted for. 

First of all, it should be noted that we can integrate Eq. over t' — t” 
because of the presence of the term e^cY' — t") in r(q, t') — r(q', t") [see Eq. 
fl29|) ]. The corresponding fast oscillating function, -t”)^ appears in the 

last factor of Eq. (jnjl . Integration of this factor results in 



t")e 


iezgc{t'—t") 



(33) 
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The lower and the upper limits of the integration over t' — t" are replaced by 
=Foo. This can be approved when the propagation time, t, is much greater 
than Lq/c. In other factors in Eq. ([3T]), the substitution t” = t' is used. 

The relation fl3^ means, that only the gx^y- components enter Eq. fl3T|) . In 
particular, the Fourier-transform ‘ip{g) should be considered as a function of 

the two-dimensional vector gj_: ip = ip { \/gl + g^ ) • This observation corre¬ 


sponds to the known Markov approximation [3] where it is assumed that the 
index-of-refraction fluctuations are delta-function correlated in the direction 
of propagation. In fact, our derivation, based on the paraxial approximation, 
supports the validity of the Markov approach which at first sight seems to 
be doubtful. 

Using Eqs. (12^ and fl33|) . the expression fiM]) is simplified to 


'jpp/ — 




f dt' [ dgiPig) g ■ (p + kt'-) g ■ (p' + 

Jo J qo' go 


(t-f) 


X (e 


/t* c«il(il-i'){F[r(q,ti)]-F[r(q',ti)]} 


). 


(34) 


where all the vectors have only the x— and ^/-components, and Cq_q/ = 

c(q- q')/go- 

As we see from Eq. fl5T|l . to obtain cpppi one needs to calculate the 
average value of the exponential function which is similar to the function in 
). Following the previous procedure, this average can be rewritten as 


exp < - tg- 


go 


dti{ti -T){F[r(q,fi)] - F[r(q',ti)]} 


= exp -27rc3 / dh{h-t'f / dgV(//')(g-g')"[l- 


(35) 

Again, the same function appears in the exponent of the right-hand side of 
Eq. (135|1 after using the trajectories (|28ll . Similar steps can be undertaken 
many times. In this way, the time hierarchy, 0 < < fi... < < f, is 

generated. If the photon-turbulence interaction time, t — ti, is short, the 
disturbance of the trajectory is small and vanishes when ti —)■ t. In this case 
both values, r(q, U) and r(q',U), approach the value of r irrespective of the 
initial momenta q and q'. Therefore we substitute the quantity 


^((g'- [r(q,ti) -r(q',fi)])^) 


(36) 
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instead of 


,-ig'-[r(q,ti)-r(q',ti)] 


) 


( 37 ) 


assuming the exponent in Eq. fl37|l to be small. The linear in g' term in 
the expansion of the exponential factor is ignored because of its zero-value 
contribution into the integral over g' in Eq. (15^ . Then the term flTr|) reduces 
to 

2 ^^ (Cq-q'-gO^ + ^(^-^l)^ / C?gX£/")(Cq-q'-g")^(g'-g")^- (38) 


To obtain Eq. (j38|l . the approximate relation, 

F„[r(q,t2))] - Fo[r(q',t2))] ~ Cq_q/(t 2 - t)drF^[r + Cq(t 2 - t)], (39) 


where ti < t 2 < t, was used. This approximation is in the spirit of the previ¬ 
ous step, where the turbulence effect was assumed as a small perturbation. 

Substitution of Eq. fl55]l into the right-hand side of Eq. fl5S]l and integra¬ 
tion over variables g', g" and ti result in 


exp 


- 2.52 ■ W~^C% 



, - ff 

560 

2-560 


(40) 


where /g = /o/27r, and 9 is the angle between the two-dimensional vectors g 
and q — q'. 

After substitution of fITOD into fl35D . fl3^ into fIMl) and fIMD into fl25|) . we 
calculate (/^(r, t)). Many-fold integrations over the variables q, q', p, p', k, k', 0, 
and t' are performed mainly numerically with employing a computer cluster. 

In the course of integration, we have used the Tatarskii modihcation of the 
refractive index spectrum which is derived from the von Karman form flT^ 
by setting Lg ^ = 0. The results for cx^ are shown in Figs. 1-3. 


4 Discussion 

Figs. 1-3 can be used to illustrate the importance of the correlations of 
different trajectories. To simplify our argumentations, we consider a coherent 


13 








laser beams, i.e., the case rg = ri. Two terms in the last braces of Eq. 
(1^ contribute equally into {P{r,t)). Moreover, if one sets (fyppi = 0 in Eq. 
fl2T|) . thus ignoring the correlations of photons with different initial momenta, 
we obtain {/^(r, f)) = 2(/(r,t))^. The scintillation index, cx^, is equal to 
unity here. This physical picture is realized for a long-distance propagation 



Figure 1: Scintillation index of a coherent and partially coherent beams in 
the atmosphere versus propagation distance 2 ;. Dashed curves correspond 
to the multiplicative approximation fl3^ for the photon correlations; solid 
curves are obtained within the present paper’s approach [see Eqs. ([35]-00])]. 

tq = 0.01m, = 10“^m, and go = 10^m“^. The upper 

two curves correspond to the coherent beam. 

it —)■ cxd) when the oscillating factor ^ conhnes the effective volume 

of the integration over g and q — q' to zero [see Eq. fl3T|) ]. For hnite values 
of t, the contribution of (j)pp' becomes quite sizeable that is seen in Figs. 1-3 
where the values of are greater than unity. There is a positive contribution 
of (j)pp' term into the last exponent in Eq. (12T|) when the vectors p, k and 
p', k' have opposite signs and the difference |q — q^| is not too large. The 
most favorable conditions are realized when 

p = -p', k = -k', q=q'. (41) 
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Figure 2: The same as in Fig. 1 but for a weaker turbulence strength: 
Cl = 2.5 X lO-i^m-2/3. 




In this case the sum (ppp + 2(pppi + (pp^p/ is equal to zero . Eqs. fHT]) can 
be interpreted as the ’’super-correlation” conditions under which the value 
of M is equal to unity and does not depend on the turbulence. 

The dependence of cx^ on the initial radius tq can be explained as follows. 
The characteristic values of the initial momentum, q ~ \/2/ro, is greater for 
small tq. Hence the volume of integration over q — q' is also greater. At the 
same time the corresponding increase of 0pp/ occurs only for short distances, 
z, where time intervals t are sufficiently small and the oscillating factor in Eq. 
fl3T|) is close to unity. Therefore, when tq decreases, there is an increase of 
accompanied with the displacement of the region with enhanced fluctuations 
towards small z. This is clearly seen in Fig. 3. 

In a similar way we can explain a considerable difference of found for 
the plane-wave and spherical-wave models of radiation in Ref. [M] (Figs. 1 
and 2 there). It follows from the above reasonings that this effect arises due 
to very different initial q-volumes in the two models. 

Also, the calculations of in the Ref. [15] should be mentioned where a 
simplified model of the turbulence was used (see Fig. 1 there). The results 
of Ref. [15] well correlate with ours. 

Comparing the results of the present paper and those, based on the ap¬ 
proximation of uncorrelated trajectories (l3^ (respectively, solid and dashed 
lines in Figs. 1 and 2), we see a more pronounced growth of at a moder¬ 
ate turbulence in the former case. Figures 1-3 illustrate that this holds true 
for the distances of 1 — 3 km. We attribute the evident distinction of the 
results to a better accuracy of accounting for the correlations of the photon 
trajectories. At the same time, both approaches provide the known in the 
literature saturation effect: —)■ rf/vg when z oo. 

The phase diffuser with a short characteristic time (a high-frequency dif¬ 
fuser) does not change qualitatively the physical picture described above. At 
the same time, both approaches reveal an ability of the diffuser to suppress 
scintillations which is favorable for communication performances. 

The effect of the phase diffuser is explained as follows. The initial phase 
relief, introduced by the diffuser, varies in time. The photon trajectories 
depend on the initial state of the radiation and varies synchronously with 
the diffuser state. A “slow” detector integrates the contribution of these 
photons. Although the atmosphere stays almost frozen during the integration 
time, the diffuser provides a better averaging of the propagating radiation 
over the refractive-index relief. Therefore, the fluctuations of the detected 
signal decrease. 
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4.0-1 



z, km 


Figure 3: Scintillation index versus propagation distance z; for different initial 
radii of the beam: ro = 0.01m, 0.03m, 0.05m. The rest of the parameters 
are the same as in Fig 1. 

This is not a unique way to suppress fluctuations. For example, the 
authors of Ref. |3S] proposed to use asymmetric optical vortices. The range 
of a weak and moderate turbulence was studied. Numerical simulations of 
the beam propagation showed promising results. It should be emphasized 
that in this case the experimental setup does not require a high-frequency 
phase diffuser. 

5 Applicability of the distribution function 
approach for short distances 

Our analysis is based on Eq. fl2ip obtained within the concept of photon 
trajectories. To consider photons as particles, whose density in the (r, q) 
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domain is defined by the distribution function /(r,q, t), the uncertainty of 
the momentum, q, should be small. The value of the uncertainty can be 
estimated from the dehnition of the distribution function ([1]) as /c/2. It 
follows from Eq. ffTTl) that close to the source and in the absence of the 

diffuser the ratio ^ ~ ^fc 2 / 4 q /2 = ( 2 /r^|i /2 = 1- Hence in the vicinity of the 
source, our calculations of are not applicable if the light is in a coherent 
state. 

The situation changes drastically for a remote detector. With increase of 
the propagation path, z, the value of q increases. The corresponding gain of 
the photon momentum, Aq, is generated by a random force, F. Hence the 
average value, (Aq), is equal to zero while the nonzero mean-square value is 
given by [12] 

(Ag2) = 0.0667r2r(l/6)go'/o-'^'C'^- (42) 

In contrast to q, the value of k decreases because of the broadening of the 
beam. The mean-square of the beam radius is given by [Slile] 

8z^T- 


R^ = i[i + 4^ 

2 q^r^rf 


+ 




( 43 ) 


where T = 0.558/g When the last term in square brackets dominates, 

the ratio (f /{k/2Y can be estimated as 
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(44) 


where (Ag^) is assumed to be of the order of g^ thus ignoring the square of 
the initial momentum 2/rQ. 

Substituting z = lO^m, gg = and Iq = 2tt ■ 10~^m into Eq. (jUj), 

we obtain ~ 21 that provides adequacy of our approach for the whole 
range of 2 ; variations shown in Figs. 1-3. This range concerns not only 
coherent, but also partially coherent beams. For partially coherent beams, 
the minimum z can be even smaller than for coherent beams. This is because 
of an additional diffuser-caused growth of g^, which is estimated by the value 

^^diffuser ~ 


6 Conclusion 

The paper continues the studies presented in Refs. [H [16]. Using the ap¬ 
proach of the distribution function, the problem of obtaining of reduces to 
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calculation of the correlations between different photon trajectories. Assum¬ 
ing the outer radius of turbulent eddies much smaller than the propagation 
distance, the iterative procedure for calculations of these correlations is de¬ 
veloped. The modihed approach makes it possible to extend applicability of 
the theory to a wider range of the propagation distances. This range includes 
a strong turbulence as well as a considerable part of a moderate turbulence 
where the scintillation index tends to reach its maximum value. The cri- 
terium, derived in Sec. 5, imposes the restriction on our theory from the side 
of short distances (weak turbulences). 
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